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Abstract 

Direct detection experiments, which are designed to detect the scattering of dark 
matter off nuclei in detectors, are a critical component in the search for the Uni- 
verse's missing matter. The count rate in these experiments should experience an 
annual modulation due to the relative motion of the Earth around the Sun. This 
modulation, not present for most known background sources, is critical for solidi- 
fying the origin of a potential signal as dark matter. In this article, we review the 
physics of annual modulation, discussing the practical formulae needed to interpret 
a modulating signal. We focus on how the modulation spectrum changes depending 
on the particle and astrophysics models for the dark matter. For standard assump- 
tions, the count rate has a cosine dependence with time, with a maximum in June 
and a minimum in December. Well-motivated generalizations of these models, how- 
ever, can affect both the phase and amplitude of the modulation. We show how a 
measurement of an annually modulating signal could teach us about the presence 
of substructure in the Galactic halo or about the interactions between dark and 
baryonic matter. We briefly review the current experimental situation for annual 
modulation, and discuss future experimental directions. 
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I. INTRODUCTION 

The Milky Way galaxy is known to be surrounded by a halo of dark matter whose 
composition remains a mystery. Only 5% of the Universe consists of ordinary atomic matter, 
while the remainder is 23% dark matter and 72% dark energy [T]. Identifying the nature 
of this dark matter is the longest outstanding problem in all of modern physics, stemming 
back to observations in 1933 by Fritz Zwicky; he proposed the existence of "Dunkle Materie" 
(German for "dark matter" ) as a source of gravitational potential to explain rapid motions 
of galaxies in the Coma Cluster [2]. Subsequently, others discovered flat rotation curves in 
disk galaxies, starting with Babcock in 1939 [3] and followed (more persuasively and with 
better data) by Rubin and Ford [I] and Roberts and Whitehurst [5] in the 1970s. Their 
results imply that the predominant constituent of mass inside galaxies must be nonluminous 
matter (see Refs. [SI Ej for reviews). 

A leading candidate for this dark matter is a Weakly Interacting Massive Particle 
(WIMP). The terminology refers to the fact that these particles undergo weak interactions 
in addition to feeling the effects of gravity, but do not participate in electromagnetic or 
strong interactions. WIMPs are electrically neutral and the average number of interactions 
with the human body is at most one per minute, even with billions passing through every 
second (8]. The expected WIMP mass ranges from 1 GeV to 10 TeV. These particles, if 
present in thermal equilibrium in the early universe, annihilate with one another so that a 
predictable number of them remain today. The relic density of these particles is 

tt x h 2 ~ (3 x 10- 26 cm 3 /sec)/(aw) ann , (1) 

where Q x is the fractional contribution of WIMPs to the energy density of the Universe. An 
annihilation cross section (av ) ann of weak interaction strength automatically gives the right 
answer, near the value measured by WMAP pQ. This coincidence is known as the 'WIMP 
miracle" and is why WIMPs are taken so seriously as dark matter candidates. Possibly the 
best WIMP candidate is motivated by supersymmetry (SUSY): the lightest neutralino in the 
Minimal Supersymmetric Standard Model (MSSM) and its extensions |9j. However, other 
WIMP candidates arise in a variety of theories beyond the Standard Model (see Refs. [10j [11] 
for a review). 

A multitude of experimental efforts are currently underway to detect WIMPs, with some 
claiming hints of detection. There is a three-pronged approach: particle accelerator, indirect 
detection (astrophysical), and direct detection experiments. The focus of this review article 
is the third option - direct detection experiments. This field began thirty years ago with the 
work of Drukier and Stodolsky [12], who proposed searching for weakly interacting particles 
(with a focus on neutrinos) by observing the nuclear recoil caused by their weak interactions 
with nuclei in detectors. Then, Goodman and Witten [13] made the important point that 
this approach could be used to search not just for neutrinos but also for WIMPs, again 
via their weak interactions with detectors. Soon after, Drukier, Freese, and Spergel [H] 
extended this work by taking into account the halo distribution of WIMPs in the Milky 
Way, as well as proposing the annual modulation that is the subject of this review. 

The basic goal of direct detection experiments is to measure the energy deposited when 
WIMPs interact with nuclei in a detector, causing those nuclei to recoil. The experiments, 
which are typically located far underground to reduce background contamination, are sen- 
sitive to WIMPs that stream through the Earth and interact with nuclei in the detector 
target. The recoiling nucleus can deposit energy in the form of ionization, heat, and/or light 



FIG. 1: A simplified view of the WIMP velocities as seen from the Sun and Earth. Due to the 
rotation of the Galactic Disk (containing the Sun) through the essentially non-rotating dark matter 
halo, the solar system experiences an effective "WIMP wind." From the perspective of the Earth, 
the wind changes throughout the year due to the Earth's orbital motion: the wind is at maximum 
speed around the beginning of June, when the Earth is moving fastest in the direction of the disk 
rotation, and at a minimum speed around the beginning of December, when the Earth is moving 
fastest in the direction opposite to the disk rotation. The Earth's orbit is inclined at ~60° relative 
to the plane of the Disk. 

that is subsequently detected. In the mid 1980s, the development of ultra-pure germanium 
detectors provided the first limits on WIMPs [To]. Since then, numerous collaborations 
worldwide have been (or will be) searching for these particles, including ANAIS [16J, ArDM 
PZj, CDEX/TEXONO [TS], CDMS [EH22], CoGeNT [23H25], COUPP [26\, CRESST [27], 
DAMA/Nal [28], DAMA/LIBRA HEED], DEAP/CLEAN Hp, DM-Ice [32], DRIFT [331 US], 
EDELWEISS [35H3JJ, EURECA [38], KIMS [39], LUX [H], NAIAD [H], PandaX 02], PI- 
CASSO g3l HE], ROSEBUD 05], SIMPLE 06], TEXONO 07J, WArP 08], XENON10 
PS], XENON100 [521 [53], XMASS [M], ZEPLIN [551 [56], and many others. 

The count rate in direct detection experiments experiences an annual modulation [hH |5"7] 
due to the motion of the Earth around the Sun (see Figure [I]). Because the relative velocity 
of the detector with respect to the WIMPs depends on the time of year, the count rate 
exhibits (in most cases) a sinusoidal dependence with time. For the simplest assumptions 
about the dark matter distribution in the halo, the flux is maximal in June and minimal 
in December. Annual modulation is a powerful signature for dark matter because most 
background signals, e.g. from radioactivity in the surroundings, are not expected to exhibit 
this kind of time dependence. The details concerning the recoil energy and modulation 
spectra depend on the specifics of both the particle physics model and the distribution of 
WIMPs in the Galaxy. We discuss these possibilities in this review. 

For more than a decade, the DAMA experiment [29J has been claiming detection of an 
annual modulation. The experiment, which consists of Nal crystals, is situated in the Gran 
Sasso Tunnel under the Apennine Mountains near Rome. By now, the amount of data col- 
lected is enormous and the statistical significance of the result is undeniable. The DAMA 
annual modulation is currently reported as almost a 9a effect |30j . and is consistent with an 
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~80 or 10 GeV [58T63J WIMP elastically scattering predominantly off of iodine or sodium, 
respectively. Many other direct detection experiments have presented null results that are 
in clear conflict with the high-mass window. The viability of the 10 GeV WIMP remains a 
controversial issue because it is not clearly compatible nor clearly incompatible with other 
null experiments once various detector systematics are taken into account. Recently, the 
CoGeNT experiment reported a 2.8a evidence for an annual modulation [25] and a third 
experiment, CRESST-II, has also announced anomalous results [27J. Whether DAMA, Co- 
GeNT, and CRESST are consistent in the low-mass window is still debated jHH [65]. Yet 
CDMS sees no annual modulation [22J, and both CDMS [20l [2TJ and XENON [SH [53] find 
null results that appear to be in conflict with the three experiments that report anomalies. 

The current experimental situation in direct detection searches is exciting. Understanding 
the anomalies and the role that different experiments play in validating them is of crucial 
importance in moving forward in the search for dark matter. In this review article, we 
seek to provide the reader with the basic tools necessary to understand a potential dark 
matter signature at a direct detection experiment, focusing on the annual modulation of 
the signal. We begin in Section [TT] by reviewing the basics of direct detection techniques for 



WIMPs, describing the particle physics in Section II A and the astrophysics in Section II B 



We describe the Standard Halo Model (SHM) as well as modifications due to substructures. 
In Section III , we examine the behavior of the annual modulation signals for both the SHM 



and substructures. In Section IV, we turn to the experimental status, briefly reviewing 
the current anomalies and null results. We conclude in Section |Vj The Appendices discuss 
quantities required for understanding results of direct detection experiments. Appendix [A] 
describes the quenching factor, and Appendix [B] presents analytical results for the mean 
inverse speed for commonly used WIMP velocity distributions, a quantity necessary for a 
computation of expected count rates in detectors. 



II. DARK MATTER DETECTION 

Direct detection experiments aim to observe the recoil of a nucleus in a collision with a 
dark matter particle (TB]. After an elastic collision with a WIMP \ of mass m x , a nucleus 
of mass M recoils with energy E nr = (fi 2 v 2 /M)(l — cos8), where fi = m x M/(m x + M) is 
the reduced mass of the WIMP-nucleus system, v is the speed of the WIMP relative to the 
nucleus, and 9 is the scattering angle in the center of mass frame. The differential recoil 
rate per unit detector mass is 

dR n Y I do \ 2p v f„ „. .da 



dE nr M \ dE nr I m x J ' dq 

where n x = p x l m x * s the number density of WIMPs, with p x the local dark matter mass 
density; /(v, t) is the time- dependent WIMP velocity distribution; and ^(q 2 ,v) is the 
velocity-dependent differential cross-section, with q 2 = 2ME nr the momentum exchange in 
the scatter. The differential rate is typically given in units of cpdkg -1 keV -1 , where cpd 
is counts per day. Using the form of the differential cross-section for the most commonly 
assumed couplings, to be discussed below, 

dR 1 



5 



where a(q) is an effective scattering cross-section and 
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(5) 



the minimum WIMP velocity that can result in a recoil energy E nT . Here, 5 is the mass 
splitting between the lightest and next-to-lightest states in the spectrum in the case of 
an inelastic scattering interaction; 1 we shall consider only the elastic scattering case for 
the remainder of this review. The benefit of writing the recoil spectrum in the form of 
Eqn. ^ is that the particle physics and astrophysics separate into two factors, a(q) and 
PxV( v mm, t), respectively. We discuss each of these factors in the following two subsections. 
More detailed reviews of the dark matter scattering process and direct detection can be 
found in Refs. HI EM] . 



A. Particle Physics: Cross-section 

For a SUSY neutralino and many other WIMP candidates, the dominant WIMP-quark 
couplings in direct detection experiments are the scalar and axial-vector couplings, which 
respectively give rise to spin-independent (SI) and spin-dependent (SD) cross-sections [9]. 
In both cases, 

^v) = ^F 2 (q)e(q m .^-q). (6) 

Here, is the Heaviside step function, g max = is the maximum momentum transfer 
in a collision at a relative velocity v, and the requirement q < g max gives rise to the lower 
limit v > f m i n in the integral for t] in Eqn. Q. In the above equation, o"o is the scattering 
cross-section in the zero-momentum-transfer limit — we shall use cxsi and <tsd to represent 
this term in the SI and SD cases, respectively — and F 2 (q) is a form factor to account for 
the finite size of the nucleus. The WIMP coherently scatters off the entire nucleus when the 
momentum transfer is small, giving F 2 (q) — > 1. However, as the de Broglie wavelength of 
the momentum transfer becomes comparable to the size of the nucleus, the WIMP becomes 
sensitive to the spatial structure of the nucleus and F 2 (q) < 1, with F 2 (q) <C 1 at higher 
momentum transfers. It is traditional to define a form-factor corrected cross-section 

a(q)=a F 2 (q), (7) 



Inelastic scattering with 5 ~ 0(100 keV) was first invoked to reconcile the DAMA anomaly with the 
CDMS limits [55] • Although this explanation has since been ruled out by XENON100 for conventional 
couplings (67], alternate formulations remain viable as a means of reconciling the experimental results 
(see e.g. Ref. [68]). More generally, inelastic scattering (for arbitrary 6) remains an interesting possibility 
for direct detection experiments, yielding distinct recoil spectra. 
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as was used in Eqn. ([3]) above. We note that this is an effective cross-section, whereas the 
actual scattering cross-section is given by f dq 2 -^(q 2 ,v) for a given relative velocity v. 

The total WIMP-nucleus scattering rate is then the sum over both the SI and SD con- 
tributions, each with its own value of the form factor. We describe these two cross-sections 
below and then briefly discuss more general operators. 



1. Spin-independent cross-section (SI) 

The SI WIMP-nucleus interaction, which occurs through operators such as (xx){qq)i has 
the cross-section 

a sl = -^\zf p + (A-Z)f a ]\ (8) 

7T L J 

where Z and A — Z are the number of protons and neutrons in the nucleus, respectively, 
and f p (/ n ) is the effective coupling to the proton (neutron). For neutralinos and most other 
WIMP candidates with a SI interaction arising through scalar couplings, f n ~ f p and the 
SI scattering cross-section of WIMPs with protons and neutrons are roughly comparable, 
Cn,si ~ Op.si- For identical couplings (/ n = / p ), the SI cross-section can be written as 



2 



trsi = ^A 2 <7 p ,si , (9) 

where fi p is the WIMP-proton reduced mass. As neutralinos are the currently favored 
WIMP candidate, this assumption is widely made throughout the direct detection literature. 
Although typically f n ~ f p , models can be constructed that violate this condition (e.g., 
isospin- violating dark matter [72\). We assume identical SI couplings for the rest of this 
review. 

The SI cross-section grows rapidly with nuclear mass. The explicit A 2 factor in Eqn. Q 
arises from the fact that the contributions to the total SI cross-section of a nucleus is a 
coherent sum over the individual protons and neutrons within. In addition, for WIMPs that 
are much heavier than the nucleus, \ ~ A 2 , so the cross-section scales as ~ A 4 . However, 

the form factor suppression becomes more significant as the size of the nucleus increases, 
so the scattering rate does not scale as ~ A A for heavy nuclei, though it still rises rapidly 
with A. As a result, direct detection experiments often use heavy nuclei to increase their 
sensitivity to WIMP scattering. 

The SI form factor is essentially a Fourier transform of the mass distribution of the 
nucleus. A reasonably accurate approximation is the Helm form factor j7H [73] : 

F(q) = 3e - 9 V/2 sin(gr n ) - gr cos(gr n ) (1Q) 

where s ~ 0.9 fm and r 2 = c 2 + ^n 2 a 2 — 5s 2 is an effective nuclear radius with a ~ 0.52 fm 
and c ~ 1.23A 1 / 3 — 0.60 fm. Further details on SI form factors can be found in Refs. [7T| TTQ. 



2. Spin- dependent cross-section (SD) 



SD scattering is due to the interaction of a WIMP with the spin of the nucleus through 
operators such as (x7^75X) (<?7 A1 75?) 5 an d takes place only in those detector isotopes with an 



unpaired proton and/or unpaired neutron. The SD WIMP-nucleus cross-section is 



32/j, 2 



G 2 F J(J +1)A 2 , 



(11) 



0"SD = 



7T 



where Gp is the Fermi constant, J is the spin of the nucleus, 




(12) 



where (S p ) and (S n ) are the average spin contributions from the proton and neutron groups, 
respectively, and a p (a n ) are the effective couplings to the proton (neutron). Unlike the 
SI case, the two SD couplings a n and a p may differ substantially (though they are often of 
similar order of magnitude), so that a simplification comparable to Eqn. ^ for SI scattering 
is not made in the SD case. Because of the uncertain theoretical relation between the two 
couplings and following from the fact that one of (S p ) or (S n ) is often much smaller than the 
other, experiments typically only significantly constrain one of the two SD cross-sections, 
a p,SD or a n sD, but not both. 

SD scattering is often of lesser significance than SI scattering in direct detection exper- 
iments for two main reasons. First, SI scattering has a coherence factor A 2 that the SD 
scattering is missing. In fact, the spin factors J, (S p ), and (S n ) are either zero or 0(1), 
so the SD cross-section does not grow as rapidly with nucleus size as the SI cross-section 
does. Thus, whereas o"si oc A A for heavy WIMPs, o"sd oc A 2 (n.b. this remaining A 2 factor 
arises from ^ ~ ^4 2 )- Second, spin- zero isotopes do not contribute to SD scattering, so 
the SD scattering is reduced in elements where non-zero-spin nuclei represent only a small 
fraction of the naturally occurring isotopes within a detector's target mass. We note that 
SD couplings may often be larger than SI couplings; e.g., for an MSSM neutralino, it is 
often the case that o" p SD / '<t Pi si ~ C(10 2 — 10 4 ). However, even with this ratio of couplings, 
SI scattering is still expected to dominate for the heavy elements used in most detectors for 
the two reasons described above. 

The SD form factor depends on the spin structure of a nucleus and is thus different 
between individual elements. Form factors for many isotopes of interest to direct detection 
experiments, as well as estimates of the spin factors (S p ) and (S n ), can be found in the 
reviews of Refs. |75j [76]. 



While scalar and axial-vector couplings are the dominant interactions for many WIMP 
candidates, such as neutralinos, they are by no means the only allowed couplings. In general, 
dark matter-nucleon interactions can be described by a non-relativistic effective theory, as 
detailed in [77J [75] . The effective theory approach is useful for highlighting the variety of 
operator interactions that can exist, and their potentially unique direct detection signatures. 

Consider momentum-dependent (MD) interactions, which can dominate for certain 
classes of theories [79H83] . In this case, the scattering rate takes the form 



3. General operators 
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where go is an arbitrary mass scale and i = SI, SD denotes whether the rate is independent 
of nuclear spin or not; is the conventional SI or SD scattering rate described previously. 
For the most commonly studied operators, (xx)(qq) and (x7/y5X)(?7 M 75<z)j n = and i = 
SI, SD, respectively. Generalizations to these scenarios include the operator (xi5X)(QQ)i 
which yields an exponent n — 1 and a rate that is not dependent on nuclear spin. In 
contrast, (X7sx) (975<7) nas n = 2 and i = SD. The momentum dependence in the rate has 
an important effect on the recoil spectrum, suppressing scattering at low energies. This 
leads to a peaked recoil spectrum and potentially more high-energy events than would be 
expected for the case of standard elastic scattering with no momentum dependence, where 
the rate falls off exponentially. 



B. Astrophysics: Dark Matter Distribution 

The velocity distribution /(v) of dark matter particles in the Galactic halo affects the 
signal in dark matter detectors. Here, we discuss the velocities of the dark matter compo- 
nents of the halo. The dominant contribution is a smooth virialized component, discussed 



in Section II B 1 The formation of the Milky Way via merger events leads to significant 



structure in both the spatial and velocity distribution of the dark matter halo, including 



dark matter streams and tidal debris, as discussed in Section II B 2 

Velocity distributions are frequently given in a frame other than the lab frame to simplify 
their analytical form. In this review, we define /(v) as the distribution in the rest frame 
of the dark matter population (i.e. the frame in which the bulk motion of the dark matter 
particles is zero); in the case of the (essentially) non-rotating smooth halo background, 
that frame is the Galactic rest frame. The lab frame distribution is obtained through a 



Galilean transformation as described in Section III More details of several commonly used 
distributions, including analytical forms for the mean inverse speed 77, can be found in 
Appendix [B} 

1. Smooth Halo Component 

The dark matter halo in the local neighborhood is most likely dominated by a smooth 
and well-mixed (virialized) component with an average density p x ~ 0.4 GeV/cm 3 . 2 The 
simplest model for this smooth component is often taken to be the Standard Halo Model 
(SHM) [IUE7], an isothermal sphere with an isotropic, Maxwellian velocity distribution and 
rms velocity dispersion a v . The SHM is written as 

{_J_ ( _jA 3/2 p" 3v2 /2^ f or |vl < 7; 
N csc {2nat) e > IOr |V| < <U esc ^ 

0, otherwise. 



2 Estimates for the local density of the smooth dark matter component are model dependent and vary in 
the literature by as much as a factor of two [84-88 1. Historically, 0.3 GeV/cm 3 has often been assumed 
when making comparisons between direct detection results. While this density is by no means ruled out 
by current observations, more recent estimates suggest a value closer to 0.4 GeV/cm 3 . Both values of the 
local density can be found in recent direct detection literature. 
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Here, 



N csc = erf (z) - 
with z = v csc /vq, is a normalization factor and 

v = \/2j3a v 



2 - z * 

ze 

7T 



(15) 



(16) 



is the most probable speed. The Maxwellian distribution is truncated at the escape velocity 
v esc to account for the fact that WIMPs with sufficiently high velocities escape the Galaxy's 
potential well and, thus, the high-velocity tail of the distribution is depleted. The dark 
matter escape velocity in the Milky Way is estimated from that of high- velocity stars. The 
RAVE survey finds that the 90% confidence range is 498-608 km/s |89J . Figure [2] shows 
the SHM speed distribution in the lab (Earth) frame, after accounting for the motion of the 
solar system relative to the Galactic rest frame, as well as the mean inverse speed r\. 

The sharp cut-off at the escape speed in Eqn. (14) is not physical. To smoothen the 



transition near the escape speed, one may use the (still ad hoc) distribution: 



/(v) 



Ac. 



271(7,2 



3/2 



-3v 2 /2<j2 



for 



V < v P 



otherwise . 



where 



erf (z) 



71 



:i+§*v 



(17) 



(18) 



In another approach, Ref. (90] uses King models to obtain the velocity distribution, handling 
the finite size and mass of the Galaxy in a more self-consistent manner. In these models, 
the probability distribution can reach zero at a lower velocity than the escape velocity; 
essentially, the highest bound velocities are unpopulated. In general, because of the large 
uncertainty in modeling the tail of the velocity distribution, one should approach any result 
that depends sensitively on high velocity predictions with caution. 

For the conventional SI and SD elastic scattering, the recoil spectrum falls off exponen- 
tially in the Galactic rest frame for the SHM (neglecting form factors), due to the exponential 
drop off with velocity in Eqn. (14). Even when form factors and the motion of the Earth 



through the halo are accounted for, the spectrum is still approximately exponential in the 
lab frame: 

.ID 

-Em/ Eq 



dR 
d~E n 



~ e 



(19) 



where Eq is some effective scale that is (9(10 keV) for typical WIMP and nuclear target 
masses, so that the largest contribution to the rate in detectors is at low recoil energies. For 
momentum-dependent interaction operators or inelastic scattering, the rate may instead 
peak at higher values of recoil energy. 



The isotropic, Maxwellian velocity distribution of Eqn. (14), intended to describe a class 



of smooth spherical halo models, is only a first approximation of the local halo profile. 
As reviewed in Ref. [HI], oblate/prolate or triaxial halos would be expected to have an 
anisotropic velocity distribution, which may be approximated as 
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where V{ are the WIMP velocities along three perpendicular directions with dispersions <jj. 
In general, changes to the halo shape from anisotropy result in 0(10%) changes in the annual 
modulation signal [92, 93J, although a more exact statement depends on the dark matter 
properties and the detector threshold. 

High resolution cosmological N-body simulations provide evidence that a Maxwellian 
distribution does not fully capture the velocity distribution of the smooth halo component, 
particularly along the high velocity tail, which is important for detection of low mass WIMPs 
as detectors are sensitive only to the highest velocity WIMPs in this case. Ref. [M] deter- 
mined the velocity distribution from two of the highest resolution numerical simulations of 
Galactic dark matter structure (Via Lactea II [951 [96] an d GHALO [97J). They found more 
low speed particles than in a Maxwellian case, and a distribution with a peak that is flatter 
in shape. Alternatively, analytic fits for producing better agreement with numerical results 
at the high speed tail have been obtained, such as the ansatz in Ref. [95] . 

Another issue is that most simulations only contain dark matter particles; simulating 
baryonic physics is extremely difficult, but important given that baryons dominate in the 
inner regions of the Milky Way. Gas cooling changes halo shapes from prolate-triaxial 
to more spherical when baryons are added [99H103J. with velocity distributions that are 
expected to deviate less from the standard Maxwellian than those found in dark matter- 
only simulations. Ref. [104] studied predictions for dark matter experiments within the 
context of an isolated numerical model of a Milky Way-like system designed to reproduce 
the basic properties of the Galaxy by including an equilibrated Galactic stellar disk and 
the associated Sagittarius galaxy impact, in addition to dark matter. The resulting dark 
matter velocity distribution still exhibits deviations from Maxwellian and the calculated 
recoil spectrum has an increased number of scattering events at large energies. 

Using cosmological simulations, Refs. |105H107| identified the possibility of a disk-like 
dark matter component ("dark disk") that forms from satellite merger events. Ling et al. 
[108J performed a high-resolution cosmological N-body simulation with baryons. They study 
a Milky Way sized object at redshift z = that includes gas, stars and dark matter to 
characterize the co-rotating dark disk, which could play an important role in direct detection 
experiments |lU9j . Equilibrated self-gravitating collisionless structures have been shown to 
exhibit Tsallis distributions [110H113| : 



/(v) 



N(v ,q) 



l-(l-g) 



v 21?/(l-9) 



(21) 



where N(vo,q) is a normalization constant and the Maxwell-Botzmann distribution is re- 
covered by taking the limit q — > 1. For a spherical shell at the same radial distance as the 
Sun in the Ling et al. simulation, the velocity distribution is best fit by a Tsallis distribution 
with vq = 267.2 km/s and q = 0.773. In an analysis of the dark matter and stars in the 
cosmological hydro simulation of Ref. |114j . the work of Ref. |115j has also found higher tan- 
gential motion in dark matter particles close to the disk plane than away from it, consistent 
with a dark disk. 



2. Unvirialized Structure of Halo 



The Milky Way halo forms through the merging of smaller dark matter subhalos. These 
merging events can lead to significant structure in both the spatial and velocity distribution 
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FIG. 2: Comparison of the Standard Halo Model (SHM) and an example stream, representative 
of the smooth background halo and a cold flow, respectively. The stream, modeled after the 
Sagittarius (Sgr) stream, is roughly orthogonal to the Galactic plane with speed ~350 km/s relative 
to the Sun. Upper Panel: The speed distribution (one dimensional f(v) in the frame of the Earth) 
for both components. Lower panel: The differential signal in a detector is directly proportional 
to the mean inverse speed rj(v m i n ). Here, the x-axis is u m ; n , the lower limit of the integration in 
Eqn. Q. The approximately exponential SHM and stepdike stream 77's are each shown at two 
periods of the year, corresponding to the times of year at which 7/ is minimized and maximized; 
note these times are different for the two components. 



of the dark matter halo. High resolution cosmological dark matter simulations, such as 
Via Lactea [951 ES], GHALO [97], and Aquarius |116j . find residual substructure from the 
merging process that includes dark matter clumps, cold streams, and debris flows. The dark 
matter affiliated with any of these substructures located in the Solar neighborhood affects 
count rates and spectra as well as the phase and amplitude of the annual modulation in 
experiments [Ml EH EIM23! • 

An example of spatially-localized substructure is a dense clump or subhalo of dark matter. 
If the Earth is sitting in such a clump, the local dark matter density would be larger 
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than currently expected, increasing scattering rates in experiments. According to numerical 
simulations, however, local density variations due to the dumpiness of the dark matter halo 
are unlikely to significantly affect the direct detection scattering rate. Based on the Aquarius 
Project, Ref. |124] reports that the dark matter density at the Sun's location differs by less 
than 15% from the average at more than 99.9% confidence and estimates a probability of 
10~ 4 for the Sun being located in a bound subhalo of any mass. The possibility that the 
Earth may reside in a local underdensity due to unvirialized subhalos throughout the Galaxy 
should also be taken into account when interpreting direct detection null results; Ref. [125J 
predicts a positively skewed density distribution with local densities as low as one tenth the 
mean value, but probably not much less than half. 

In addition to structure in configuration space, the dark matter halo can also exhibit 
velocity substructure in the form of debris flows or cold tidal streams. Debris flows are 
an example of spatially-homogenous velocity substructure that consists of the overlapping 
shells, sheets, and plumes formed from the tidal debris of the (sum total of) subhalos falling 
into the Milky Way |123[ I126j . Although this dark matter component is spatially uniform, 
the distribution of its Galactocentric speeds is roughly a delta function. 3 In Via Lactea II, 
more than half of the dark matter near the Sun with (Earth-frame) speeds greater than 
450 km/s is debris flow. At higher speeds, debris flow comprises over 80% of the dark 
matter. As a result, debris flow is particularly important for experiments that probe the 
high velocity tail of the dark matter distribution, such as searches for light dark matter or 
experiments with directional sensitivity. 

Tidal streams are another unvirialized component of the halo and also consist of material 
stripped from infalling satellites. As the material in the stream has not had the time to spa- 
tially mix, the stream has a small velocity dispersion in comparison to that of the virialized 
halo. A dark matter stream is coherent in velocity space, with 

/str(v) = 5 3 (v) (22) 

in the limit of zero dispersion. In some cases, particularly when examining the annual 
modulation signal, it may be important to account for the small but non-negligible dispersion 
of the stream. In such Maxwellian velocity distribution can be used, albeit with 

a much smaller a v than that for the SHM. The speed distribution of an example stream in 
the lab frame, as well as the corresponding rj, are shown in Figure [2] 

The Sagittarius stream is one of the most stunning examples of a stellar stream in our 
Galaxy. Sagittarius (Sgr) is a satellite galaxy that is located inside the Milky Way on 
the opposite side of the Galactic Center from the Sun; it is currently being disrupted and 
absorbed by the Milky Way. The Sloan Digital Sky Survey and the Two Micron All Sky 
Survey [127-129J have traced the stellar component of the tidal stream |130[ 1131 j . Two 
streams of matter are being tidally pulled away from the main body of the Sgr galaxy 
and extend outward from it. Whether the Sgr stream passes close enough to the solar 



3 Note the distinction between a delta function in the speed for a debris flow and a delta function in the 
velocity for a stream (below). 

4 Tidal streams can have much more anisotropic velocity distributions than the smooth halo background, 
with a larger dispersion along the longitudinal direction than transverse directions 118 j. Still, the isotropic 
Maxwellian distribution with an appropriately chosen o~ v can provide a sufficiently good approximation 
for the purposes of direct detection calculations. 
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neighborhood to affect direct detection experiments remains up for debate. Early data 
indicated that the leading tail of stellar material ripped from the Sgr galaxy passes only a 
few kpc from the solar neighborhood |132[ 1133] . but later studies indicated that the center 
of the stream's stellar component could be farther away than initial estimates suggested 
|134j . Most recently, however, Ref. |1U4] analyzed self-consistent N-body simulations of the 
Milky Way disk and the ongoing disruption of the Sgr dwarf galaxy and argued that the 
dark matter part of the Sgr stream may, in fact, impact the Earth. Streams can have a 
variety of effects on direct detection experiments \119\ I120j . as discussed further below. 

Alternative models of halo formation, such as the late-infall model [135H137] more recently 
examined by Sikivie and others [138H144] . also predict cold flows of dark matter. In the 
caustic ring model |145j . the annual modulation is 180 degrees out of phase compared to the 
usual (isothermal) model. Any such streaming of WIMPs (we will henceforth use "stream" 
to imply any cold flow) will yield a significantly different modulation effect than that due to 
a smooth halo. 

Finally, there may be unbound dark matter of extragalactic origin passing through the 
Galaxy. If present, these high-speed particles can increase the number of high-energy scat- 
tering events in a direct detection experiment [146L I147j . 



III. ANNUAL MODULATION 

The velocity distribution in the Earth's frame /(v,t) changes throughout the year due to 
the time- varying motion of an observer on Earth. Assuming /(v) is the velocity distribution 
in the rest frame of the dark matter population, i.e. the frame where the bulk motion is 
zero, the velocity distribution in the lab frame is obtained after a Galilean boost: 

/(v,t) = 7(v obB (t)+v), (23) 

where 

VobsOO = v +V ffi (t) (24) 

is the motion of the lab frame relative to the rest frame of the dark matter, v is the 
motion of the Sun relative to that frame, and V (t) is the velocity of the Earth relative 
to the Sun. For a non-rotating, smooth background halo component, such as the SHM, 
v = vlsr + v 0i p ec , where vlsr = (0,t> ro t,0) is the motion of the Local Standard of Rest 
in Galactic coordinates, 5 and v pec = (11,12,7) km/s is the Sun's peculiar velocity (see 
e.g. Refs. |148[ 1149] and references therein). The canonical value for the disk rotation speed 
v rot has long been 220 km/s |150j . but more recent estimates tend to place it 5-15% higher 
[T5THT53] . 



The V (t) term in Eqn. (24) varies throughout the year as the Earth orbits the Sun, 
leading to an annual modulation in the velocity distribution and, thus, the recoil rate. 
Written out in full, 

V©(t) = V® [ii cosuj(t — ti) + i 2 sinu;(i — ti)] , (25) 



5 Galactic coordinates are aligned such that x is the direction to the Galactic center, y is the direction of 
the local disk rotation, and z is orthogonal to the plane of the disk. 
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where u = 27r/year, V© = 29.8 km/s is the Earth's orbital speed around the Sun, and £\ 
and £2 are the directions of the Earth's velocity at times ti and t\ + 0.25 years, respectively. 



Equation (25) neglects the ellipticity of the Earth's orbit, which is small and only gives neg- 
ligible changes to the velocity expression (see Refs. [7TI 1154] for more detailed expressions). 
In Galactic coordinates, 

e x = (0.9931, 0.1170, -0.01032) and e 2 = (-0.0670, 0.4927, -0.8676) , (26) 

where E\ and £2 are the directions of the Earth's motion at the Spring equinox (March 21, 
or ti) and Summer solstice (June 21), respectively. 

If we define the characteristic time to as the time of year at which v b s (t) is maximized, 
i.e. the time of year at which Earth is moving fastest with respect to the rest frame of the 
dark matter, then the magnitude of v Q b s (t) is approximately 

Vobsit) ^v G + bV® cos u (t -t ), (27) 

where b = ^/bf + b\ for b { = Si ■ v is a geometrical factor associated with the direction 
of v bs relative to Earth's orbital plane (note |6| < 1). The approximation is valid when 
V®/vo -C 1, as is the case with nearly all halo components. 



Once the Galilean transformation of Eqn. (23) is performed, r](v mm ,t) is calculated via 
Eqn. Q; see Appendix [B] for analytical forms of 7] for several commonly used distributions. 
Because the modulation rate must have a fixed period of one year, 6 the differential scattering 
rate can be expanded in a Fourier series: 

dR °° 00 

^Omin, t) = + A n cosnu(t - t ) + ^ B n sinnw(t - t ) , (28) 

n=l n=l 

where the Fourier coefficients A n and B n are functions of f m in- If the velocity distribution 
in the rest frame of the dark matter is isotropic, then B n = 0. This simplification is a direct 
result of expanding about t : though the Fourier expansion could be made about any other 
(arbitrary) phase, the sum would include both cosine and sine terms. The expansion in 
terms of only cosines should not be surprising as f bs(£) contains only a single cosine term 



in Eqn. (27). This simplification does not apply for anisotropic distributions; however, for 
nearly all realistic anisotropic distributions, B n A n , and the sine terms in the expansion 
can still be neglected. As a consequence, the modulation will always be symmetric (or very 
nearly so) about the characteristic time to- 

For smooth components of the halo, such that f(v) is slowly varying over <5v ~ Vq, one 
further finds that A ^> A\ ^> A n > 2 , assuming v Q ^> V e as is the case for most components 
of the halo. This relation further holds for any structure in the halo when /(v) is slowly 
varying for |v| ~ v m [ n , as is the case for cold flows where v m m ^ v & . Higher order terms 
in the Fourier expansion may become important when /(v) exhibits sharp changes in the 
vicinity of |v| ~ f m i n , which can happen in the case of a stream. 



6 The density and intrinsic velocity distribution will change as the solar system moves into, through, and 
back out of any substructure of hnite size, such as a clump, leading to variations in the recoil rate that 
do not manifest as an annual modulation. However, the time scales involved are typically many orders of 
magnitude longer than a year and such temporal variations can be neglected. 
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SHM Debris Flow Stream 




FIG. 3: A comparison of the shapes of the total rate shown at two periods of the year, correspond- 
ing to the times of year at which the rate is minimized and maximized, as well as the modulation 
amplitude, for three different halo components: SHM (left), debris flow (middle), stream (right). 
The normalization between panels is arbitrary. 

Except for the special cases described above, the annually modulating recoil rate can be 
approximated by 

—(E,t)nS (E)+S m (E)cosuj(t-to), (29) 

with \S m \ <C So, where S is the time-averaged rate, S m is referred to as the modulation 
amplitude (which may, in fact, be negative), u = 27r/year and t is the phase of the mod- 
ulation. 7 The quantities So and S m correspond to A and A\, respectively, in the Fourier 



expansion of Eqn. (28), but the former are the standard notation in the literature when only 
the constant and first cosine term of the Fourier expansion are used. 

In the remainder of this section, we examine the modulation for the SHM and substruc- 
ture components. Figure [3] summarizes the conclusions we reach. Note that the expected 
modulation amplitude depends sensitively on the assumed dark matter velocity distribution. 
In reality, the local dark matter is likely comprised of both a virialized and unvirialized 



7 Experiments may quote the average amplitude over some interval, 
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component, meaning that a signal at a direct detection experiment may be due to several 
different dark matter components. In this case, a modulation of the form given by Eqn. (29) 
with a fixed phase to may not be a good approximation; the shape of the modulation for 
the total rate may no longer be sinusoidal in shape and/or the phase may vary with i> m i n . 
Furthermore, there are cases when Eqn. (29) is a bad approximation even for a single halo 
component; an example will be shown below for a stream. We close this section with a 
discussion of what can be learned about the local halo in these more complicated scenarios. 



A. Smooth Background Halo: Isothermal (Standard) Halo Model 

We now apply our general discussion of the modulation rate to the example of a simple 
isothermal sphere [57]. As discussed in the previous section, the SHM is almost certainly 
not an accurate model for the dark matter velocity distribution in the Milky Way. However, 
its simple analytic form provides a useful starting point for gaining intuition about the 
modulation spectrum of the virialized dark matter component. 

As we showed in Eqn. (|3]), the differential count rate in a detector is directly proportional 
to the mean inverse speed rj; the time-dependence of the recoil rate arises entirely through 
this term. To study the expected time-dependence of the signal in the detector, we therefore 
focus on the time-dependence of rj; in particular, we can investigate the annual modulation 
of the quantity rj as it is the same as that of the dark matter count rate. 

For the SHM or any dark matter component with a velocity distribution described by 



Eqn. (14) or Eqn. (17), the mean inverse speed has an analytical form, presented in Ap- 
pendix]B| and in Refs. |121[ I155j . Figure [2] illustrates 77(f m in) for the SHM, taking t> = v TOt 
as expected for an isothermal spherical halo. 

Figure [2] shows r](v m - m ) at t ~ June 1, the time of year at which Earth is moving fastest 
through the SHM, as well as on December 1, when the Earth is moving slowest; there is a 
(small) change in r\ over the year. The corresponding recoil spectra, as a function of recoil 
energy, are given in schematic form in the first panel of Figure [3] The amplitude of the 
modulation, 



ME) ~ I 



dR , _ , . dR , _ _ . 
-(£, June 1) - _(£, Dec 1) 



(31) 



is also shown in the figure. Two features of the modulation are apparent for the SHM: (1) 
the amplitude of the modulation is small relative to the average rate, with an exception to 
be discussed below, and (2) the amplitude of the modulation changes sign at small v m i n (low 
recoil energies). This phase reversal can be used to constrain the WIMP mass. 

Figure [4] illustrates the residual time- varying signal for the SHM (dashed blue curves). 
The different panels show how the modulation depends on v m \ n . In general, the modulation 
has a sinusoidal shape and is symmetric about t ; the sinusoidal shape allows for the use 



of the amplitude approximation given by Eqn. (31). For small t> m in (low recoil energies), 
the rate is minimized at a time to, while for large t> m in (high recoil energies), the rate is 
maximized at to. 

An important quantity of interest is the modulation fraction, defined as the size of the 
modulation amplitude relative to the average total rate: S m /So. For a wide range of f m i n , 
the modulation fraction is (9(1 — 10%), as seen in Figure |lj For t> min above ~ 500 km/s, 
both So and S m fall rapidly with t> m in as scatters come only from WIMPs in the tail of the 
Maxwellian distribution. In this region, Sq falls more rapidly, so the modulation fraction 
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grows, going from 0(10%) to (9(100%). Because of the low absolute rates at these higher 
energies, experiments are generally not sensitive to the v m in region with high modulation 
fraction. However, for WIMPs that are much lighter than the nuclear target, large modu- 
lation fractions do correspond to the recoil energies of interest in detectors, and an order 
unity modulation can be observed. 8 

As noted above, the phase reversal of the annual modulation, illustrated in Figure [4], 
can be used to determine the WIMP mass [156] and is perhaps the best feature of a direct 
detection signal for doing so. While the phase of the modulation is fixed for a given f m i n , 
regardless of the WIMP mass, the phase of the modulation for a given recoil energy E m 
is not, as the E nr associated with a given t> min is dependent on the WIMP mass through 
Eqn. (|5]). Thus, an experimental determination of the recoil energy at which the phase 
reverses, which occurs at v m \ n ~ 210 km/s for the SHM, can constrain the WIMP mass. 
For a germanium detector, the SHM phase reversal is expected to occur at recoil energies 
of 1, 5, and 15 keV for WIMP masses of 10, 25, and 60 GeV, respectively; the modulation 
spectrum should be readily distinguishable between these cases. As the WIMP gets much 
heavier than the target nucleus, the recoil spectrum becomes degenerate and the energy of 
the phase reversal approaches a fixed value (62 keV for germanium); observation of a reversal 
at this energy allows only a lower limit to be placed on the WIMP mass. 9 We emphasize 
that detection of this phase reversal could constitute an important signature of a WIMP 
flux, as backgrounds would not give rise to such an effect. 

As discussed in the previous section, the SHM may not be an accurate model of the 
smooth background halo. However, the generic features of the SHM modulation signal 
discussed here are also features to be expected of any smooth background halo model. In 
particular, the modulation should be sinusoidal in shape with a phase around the beginning 
of June, have an Oil — 10%) modulation amplitude (except at high t> m in), and have a phase 
reversal [92| l93| I159j . However, all these features can be significantly altered if there is any 
significant substructure present in the halo, so we turn to substructure next. 



B. Halo Substructure 



Next, we consider the modulation spectrum when the dark matter scattering in the 
detector is dominated by unvirialized components in the halo, such as debris flows or streams. 
The schematic in Figure [3] illustrates how the modulated and total rates in these cases can 
differ drastically from that expected from a smooth halo contribution. For a complete 
discussion of the modulation spectra for dark matter streams, see Ref. |121j . A simple 
analytic approximation for the mean inverse speed of the debris flow is given in |123j . from 



8 When the size of the variations in the recoil rate throughout the year becomes comparable to the average 
rate, i.e. the relative modulation amplitude is of order unity, Eqn. (31 ) is no longer a good approximation 
to the modulation. A significant deviation from a cosinusoidal modulation would be an expected signature 
for large modulation fractions. 

9 A caveat regarding extracting limits on the WIMP mass from the phase reversal is in order. As we will 

illustrate shortly, cold flows (such as streams or caustics) can strongly affect the phase of the modulation. 
Thus, the phase of the annual modulation constrains the WIMP mass only when the distribution of 
particle velocities in the solar neighborhood is known. One may however use the results of two different 
experiments to constrain the mass without assuming a form for the velocity distribution [1571 1158) . 
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FIG. 4: The residual rate for the SHM (dashed) and an example stream (dotted) is plotted 
at several values of u m ; n . The stream is modeled after the Sgr stream. Also shown is the total 
SHM+stream modulation, assuming the local density of the stream is 10% that of the SHM. The 
residual rates are given relative to the average rate in each case, i.e. curves show the fractional 
modulation, except for the stream, where the relative rate has been divided by ten for visual clarity 
as its relative modulation is much larger. The corresponding recoil energy for each v m i n , given by 
Eqn. (jSJ), depends on the WIMP and nuclear target masses; for a WIMP mass of 60 GeV and a 
germanium target, the v m \ n values of 150, 350, and 550 km/s correspond to recoil energies of 7, 
40, and 100 keV, respectively. The phase reversal of the SHM component, which occurs at smaller 
Vxom-i can be seen by comparing the dashed blue curves in the top two panels. The recoil energy at 
which this phase reversal takes place can be used to determine the WIMP mass |156j . 
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which it is straightforward to derive the modulated and total rates (see also Appendix [B| . 

Figure [3] emphasizes that substructure components can increase the number of expected 
scattering events at large recoil energies, in comparison to the smooth halo contribution. 
This is due to the fact that velocity substructure is most likely to be found near the escape 
velocity, where the dark matter is predominantly not in equilibrium. In particular, while 
the modulation amplitude for a smooth halo falls with energy, that for debris flows and 
streams can peak at large values. Therefore, a larger-than-expected modulation at high 
recoil energies can be an important indicator for dark matter substructure in the local 
neighborhood. 

As an example of a substructure component, we now consider the case of a dark matter 
stream more fully. For the case of a dispersionless dark matter stream, the recoil spectrum 
is proportional to 

^str( u min) t) = 77\ > (32) 

where 9 is the Heaviside function, and is flat up to the cutoff energy 

E c {t) = 2 -^v ohs {tf ■ (33) 

This characteristic energy is the maximum recoil energy that can be imparted to the nucleus, 
and is obtained as follows: The maximum momentum transferred from a WIMP to a nucleus 
occurs when the nucleus recoils straight back and is q max = 2 fiv \> s (t). The maximum recoil 
energy of the nucleus then follows as E c (t) = g^ ax /(2M). A small, but non-zero velocity 
dispersion a v expected in e.g. tidal streams can soften the sharp edge of the step-shaped rj. 
The velocity dispersion for the Sagittarius stellar stream, for example, is roughly O(20 km/s) 
[127[ I129[ I131[ [160J . The dark matter in a tidal stream can be expected to have a velocity 
dispersion of a similar magnitude, though how closely it matches the dispersion of the stars 
remains an open question. 

We take as an example a stream with velocity and direction similar to what may be 
expected if the Sgr stream is accompanied by a broader stream of dark matter that passes 
near the solar neighborhood |104[ 1119^ H2Uj . This is intended as a concrete example of a 
more general phenomena, and will illustrate the basic features of an annual modulation 
signal in the presence of a stream. The stream in our example is roughly orthogonal to the 
Galactic Plane and moves at a speed ~350 km/s relative to the Sun. Its speed distribution 
and mean inverse speed r\ are shown in Figure [2] in the latter case, the step-like spectrum 
(with a softened edge) is evident. For this stream, w bs(^) is maximal on December 29 and 
minimal on June 30. 

There are two distinct features of a stream's recoil spectrum that modulate: (1) the 
height of the step and (2) the location of the step edge. Unlike the SHM, the relative 
modulation amplitude is fairly uniform at all v m[n below the velocity of the stream in the 
lab frame, though, like the SHM, the modulation is small compared to the total rate. This 
modulation, seen in the top panel of Figure |4j is sinusoidal and peaks in late December. 
Above the stream's velocity in the lab frame, the signal vanishes because there are no 
available dark matter particles (lower panel of Figure [1]). 

The modulation becomes interesting when the minimum scattering velocity is approxi- 
mately equal to the stream's velocity in the lab frame, v m i n ~ % ~ 350 km/s. This occurs 
near the edge of the step in the recoil spectrum. As illustrated in Figure [3j rj changes by 
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a relatively large amount throughout the year for t> min m v®, leading to a very large mod- 
ulation. The modulation at the softened edge of the example stream, shown in the middle 
panel of Figure |4j has a relative amplitude of nearly 70%, far larger than that for the SHM. 
Though not apparent in this figure, this is one of the special cases where the higher order 
terms of Eqn. (28) can become important and the modulation deviates from a sinusoidal 
shape; see Figure 7 of Ref. |121j for a clearer example. This very large modulation, which 
occurs only over a narrow range of recoil energies, has a phase reversed from that at lower 
energies. 

The features of the modulation for our example stream are expected of any cold flow 
(stream). Up to some cutoff energy, the modulation is uniform, relatively small, and si- 
nusoidal. Above the cutoff energy, the modulation, as well as the total rate, is negligible. 
Over a narrow range of energies about the cutoff energy, the modulation can be very large 
and possibly non-sinusoidal. However, the phase of the modulation and cutoff energy can 
vary significantly depending on the direction and speed of the stream. Observation of un- 
expected phases in the modulation and/or a narrow energy range containing an unusual 
modulation behavior would not only indicate that a significant stream of dark matter is 
passing through the local neighborhood, but would allow the direction and/or speed of that 
stream to be determined. However, there may be more than one significant stream or other 
substructure — in addition to the smooth halo background — so more than one halo com- 
ponent may make significant contributions to the recoil spectrum and modulation signals, 
complicating the interpretation of modulation data. We turn to a multiple component halo 
next. 



C. Multiple Component Halo 

Thus far, we have considered the modulation spectra of individual components of the 
dark matter halo separately. However, in greater likelihood, the local dark matter will 
be comprised of a combination of virialized and unvirialized components. In such cases, 
low-velocity dark matter will most likely be in equilibrium and well-described by a smooth 
halo, while the high-velocity tail of the distribution may have additional contributions from 
streams or debris flows. The resulting modulation spectrum will be a linear combination 
of the spectra shown in Figure [3j appropriately weighted by the relative density of each 
component. A general study of dark matter detection in the presence of arbitrary streams or 
debris flow, in combination with a smooth halo background, can be found in Refs. [121^ 123] . 

The addition of even a small amount of substructure to the smooth halo background can 
significantly alter the observational signals in direct detection experiments from those due 
to the background distribution alone. We take, for example, the case of the SHM with the 
addition of our example (Sgr- motivated) stream at a density 10% that of the SHM. While the 
overall recoil spectrum is approximately exponentially falling due to the large contribution 
from the SHM, a noticeable drop-off in the spectrum appears around a characteristic energy 
E c corresponding to the step in the stream contribution. The impact on the modulation, 
shown in Figure |4| is even more pronounced. As can be seen in the middle panel, the shape 
can differ significantly from that due to either contribution alone: the modulation is no 
longer sinusoidal and is not even symmetric in time. The phase also differs significantly: 
the peak of the modulation occurs at a time several months different from that of either 
component. From the different panels, it is also clear that the phase changes with recoil 
energy by more than just a 180° phase flip. 
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In the general case where multiple components contribute significantly to the scattering, 
the following features in the modulation spectrum can arise: 

• The phase of modulation can vary strongly with recoil energy and not just by a 180° 
phase reversal; 

• The combined modulation may not be sinusoidal, even if the modulation of each 
individual component is; 

• The combined modulation may not be time-symmetric, even if the modulation of each 
individual component is; and 

• The minimum and maximum recoil rates do not necessarily occur 0.5 years apart. 
More quantitatively, the time-dependence of the rate is no longer dominated by the A\ term 



in the Fourier expansion of Eqn. (28), and other terms in the expansion contribute. A power 
spectrum of the modulation is very useful for understanding the relative strengths of these 
higher-order contributions (see e.g. Ref. |161j ). The DAMA experiment is currently the only 
one with enough data to have produced a power spectrum of their results; their measured 
limit on A2/A1 can already provide constraints on certain types of streams, as has been 
shown for the case of inelastic dark matter in Ref. |122j . 



IV. EXPERIMENTAL STATUS OF ANNUAL MODULATION 

In this section, we discuss the current experimental status of dark matter annual mod- 
ulation searches. An extremely diverse set of direct detection experiments exists, which 
take advantage of a variety of target materials and background rejection techniques. The 
advantage of such diversity is that different targets are more or less sensitive to different 
types of dark matter and/or features in the velocity profile. For example, searches for spin- 
dependent interactions require the use of targets with non-zero spin. Also, a lighter target, 
such as germanium or sodium versus xenon, is better for detecting light mass dark matter. 

The current anomalies from DAMA, CoGeNT, and CRESST have engendered a great 
deal of excitement in the field, with debates as to whether they represent the first direct 
observation of dark matter. We now review these experiments, as well as their counterparts 
that report the tightest constraints. 



A. Current Experiments and Results 

1. The DAMA Experiment 

The Italian Dark Matter Experiment (DAMA) consists of 250 kg of radio pure Nal(Tl) 
scintillator. DAMA/Nal [2H] was the first experiment to claim a positive dark matter signal; 
it was later replaced by DAMA/LIBRA [2HJ ED], which confirmed the results. The exper- 
iment has now accumulated 1.17 ton-yr of data over 13 years of operation and claims an 
8.9cx annual modulation with a maximum in late May/early June, consistent with the dark 
matter expectation (see Figure [5j). The modulation amplitude from 2-10 keVee, taken from 
Ref. [30] , is reproduced in the top panel of Figure [6] The horizontal axis is given in terms 
of the electron-equivalent energy E ee in units of keVee, which is related to the nuclear recoil 
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FIG. 5: The residual rate measured by DAMA/Nal (red circles, 0.29 ton-yr exposure over 1995- 
2002) and DAMA/LIBRA (blue triangles, 0.87 ton-yr exposure over 2003-2010) in the 2-6 keVee 
energy interval, as a function of time. Data is taken from Refs. [2<$\ 130). The solid black line is the 
best fit sinusoidal modulation A cos[^- (t — to)] with an amplitude A = 0.0116±0.0013 cpd/kg/keV, 
a phase to = 0.400 ±0.019 yr (May 26 ± 7 days), and a period T = 0.999 ±0.002 yr [30]. The data 
are consistent with the SHM expected phase of June 1. 



energy E nr by a multiplicative quenching factor as discussed in Appendix |XJ The modula- 
tion amplitude cannot be given in terms of the nuclear recoil energy in a model-independent 
way because the experiment does not distinguish between sodium and iodine recoils on an 
event-by-event basis and the recoil energy corresponding to a given electron-equivalent en- 
ergy, related by the nucleus-dependent quenching factor, differs between the two nuclei. In 
Figure [6j the presence of a modulation is apparent below ~6 keVee, while the data above 
~6 keVee are consistent with zero modulation amplitude. 

Two possible WIMP masses can reasonably reproduce the observed modulation amplitude 
spectrum: m x ~ 10 GeV (where sodium recoils dominate) and m x ~ 80 GeV (where iodine 
recoils dominate) [58ti63] . The predicted modulation spectrum for the two best-fit masses 
and cross-sections, assuming the SHM and Si-only scattering, are shown in Figure |6j The 
behavior of the amplitude below the current 2 keVee threshold differs for the two WIMP 
masses, with the amplitude of the heavy candidate going negative. This is the phase reversal 
feature in the annual modulation discussed in Section IIII A[ Future iterations of the DAMA 
experiment, which are expected to have lower thresholds, should be able to distinguish these 
two possibilities. 

The results of the DAMA experiment are in apparent contradiction with the null results 
from other experiments, as will be discussed below. Other conventional explanations for 
DAMA's observed annual modulation have also been proposed |162j . including radon con- 
tamination and neutrons [163J. The modulating muon flux has been studied as a potential 
contaminant in the experiment |161[ I162[ I164H166] . Thus far, most of these explanations 
have been discounted (see Ref. |167j for a refutation of muons as a significant contaminant), 
but uncertainty remains. 
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FIG. 6: The bin-averaged modulation amplitude observed by DAMA (top), CoGeNT (bottom, 
red), and CDMS (bottom, blue) as a function of energy. Boxes indicate the la uncertainty for 
each bin. The DAMA results are fit at the SHM expected phase with a peak on June 1, while the 
CoGeNT and CDMS bins are given at the CoGeNT best-fit phase with a peak on April 16. The 
DAMA data are from Ref. |30| while the CoGeNT and CDMS binning and results are taken from 
Ref. |22j. To allow for direct comparison between the two germanium experiments, we present 
both CoGeNT and CDMS data in keV (nuclear recoil energy); on the other hand, DAMA data 
are presented in keVee (electron-equivalent energy) . Also shown for DAMA are the best-fit spectra 
to the data for spin- independent (SI) scattering, corresponding to a WIMP with mass 11 GeV 
(76 GeV) and SI cross-section 2xl(T 4 pb (1.5xl(T 5 pb). 

2. The CoGeNT Experiment 

The CoGeNT experiment, located in the Soudan mine in Minnesota, consists of 440 grams 
of PPC Germanium detectors with a 0.4 keVee energy threshold that makes it particularly 
well-suited to look for light dark matter [23]. Based upon 56 days of exposure, the col- 
laboration reported an excess of low energy events above the well-known cosmogenic back- 
grounds [21], which could be consistent with a ~10 GeV WIMP [1681 1169] . After more than 
a year of data-taking, an annual modulation was reported at 2.8<r with a best-fit phase of 
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April 16 [25] (see Ref. |170] for a Bayesian analysis). The lower panel of Figure [6] shows 
the modulation amplitude observed in the CoGeNT experiment for several energy bins, as- 
suming the best-fit phase; energies have been converted from electron-equivalent to nuclear 
recoil enegies as described in Appendix [Aj A significant modulation is present above 5 keV, 
which is incompatible with the total rate measured below 4 keV for standard assumptions 
about the halo and scattering properties [65] . However, the modulation could be explained 
by local substructure [&U [65], II 71 j . 

3. The CDMS Experiment 

The CDMS experiment also consists of germanium and is located in the Soudan mine. 
Using the ratio of two signals observed in an interaction with the detector target — phonons 
and ionization — CDMS can distinguish between nuclear recoil events (WIMP and/or neutron 
interactions) and electron recoil events (beta and gamma interactions), where the latter 
represents an otherwise dominant background contribution [19]. The conventional low- 
background analyses in CDMS, the most recent having 612 kg-days of exposure [20], have 
failed to detect any excess events inconsistent with background. 

To improve sensitivity to light WIMPs, which produce only low energy recoils, CDMS 
has also performed a low-energy analysis [21J, reducing its threshold from 10 keV to 2 keV. 
At these lower energies, it is more difficult to discriminate between potential signal events 
and background events, so there is far more background contamination in this analysis than 
the conventional case. The ~500 low-energy events found in this analysis are consistent 
with rough background estimates. However, CDMS places conservative no-background- 
subtraction constraints on light WIMPs, neglecting background contributions and allowing 
any or all of the events to be due to WIMPs. CDMS has also performed a modulation 
search in their low-energy data [22], finding no evidence for modulation. Constraints on the 
modulation amplitude, assuming the CoGeNT best-fit phase of April 16, are shown in the 
lower panel of Figure [6] A direct comparison can be made with the CoGeNT modulation 
results because both experiments have a germanium target, although the CDMS modulation 
analysis was only performed down to 5 keV, whereas the CoGeNT modulation data goes to 
a much lower ~2 keV. 

4- The CRESST Experiment 

The CRESST experiment, developed at the Max Planck Institute in Munich and deployed 
in the Gran Sasso Tunnel, has 730 kg-days of data with a CaW04 scintillating crystal target, 
and measures both light and heat to reject electron recoils. It reports an excess of low energy 
events with a statistical significance of over 4<r [21]. The experiment is not background-free 
however, and has experienced problems with energetic alpha and lead ions produced in 
the decay of polonium, itself produced from radon decay. For various technical reasons 
not discussed here, polonium deposited in the clamps holding the detectors in place is the 
major source of such backgrounds. The expected number of these events, which occur 
on the surface, is determined by extrapolating from high energy observations (where such 
background events are readily identifiable) to the signal regions at lower energies using Monte 
Carlo simulations. Questions have been raised as to whether the Monte Carlo simulations 
underestimate the background contamination by failing to account for the roughness of the 
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surface at microscopic scales |172j . An upcoming redesign should eliminate this background 
source, and future CRESST runs should clarify the origin of the current excess. 



5. The XENON Experiment 

The XENON collaboration has developed a series of liquid xenon target experiments, 
with the most recent iteration (XENON100) containing ~100 kg of xenon [52]. As with 
CDMS and CRESST, XENON uses two signals — scintillation and ionization in this case — 
to discriminate between nuclear recoils and electron recoils. Similarly to CDMS, XENON 
has performed both a conventional low-background analysis |53J and a low-energy analysis 
|51j . The latter sacrifices background discrimination to improve sensitivity to light WIMP 
masses. Neither analysis found an excess of events above expected backgrounds and the low- 
background analysis currently places the most stringent constraints on the SI cross-section 
for WIMPs heavier than ~10 GeV. 



6. Other Experiments 

The experiments we discussed in detail here represent only a fraction of the current direct 
detection program. No other experiment claims an excess of events consistent with dark 
matter and, for standard assumptions, none provide constraints as stringent as those from 
CDMS and XENON100. One exception is the case of SD scattering where the coupling to 
the neutron is suppressed relative to the proton (a n <C a p ). The proton-even target materials 
in CDMS and XENON only couple weakly to the WIMP in this case (see Eqn. ( 11 )), so these 



experiments place relatively weak constraints. For this case, SIMPLE [46J and COUPP 
provide the best limits. 



B. Compatibility of Experimental Results 

Figure [7] summarizes the current status of anomalies and limits for SI scattering, assum- 
ing the SHM with vq = 220 km/s and t> esc = 550 km/s. Note that the experimental limits 
and anomalies, as shown in this figure, are highly dependent on the assumptions made about 
the particle and astrophysics (see |174H177j for astrophysics- independent comparisons). The 
compatibility may change for e.g. a different WIMP-nucleus effective operator or for addi- 
tional substructure contributions. However, note that changes to the particle physics and/or 
astrophysics may change the interpretation of individual results without actually affecting 
the compatibility among different results. 10 In addition, various systematic issues regarding 
the behavior of individual detectors, such as the calibration of the recoil energy scale, can 



Newer measurements of the Sun's velocity relative to the Galactic halo (as high as 250 km/s [151] . as 
opposed to the canonical 220 km/s in common use) shifts the best fit regions and the limit curves to 
the left. For the SHM, the regions compatible with DAMA, CoGeNT, and CRESST move down in mass 
by a few GeV [178] . Because the bounds from the null experiments move to lower masses as well, the 
discrepancy between experiments is not alleviated. 
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FIG. 7: WIMP mass and SI cross-sections consistent with the anomalies seen by DAMA, CoGeNT, 
and CRESST, as well as constraints placed by the null results of CDMS and XENON. The halo 
model is assumed to be the SHM with the given parameters. The lack of overlap between the regions 
of the three anomalous results and their locations above the exclusion curves of CDMS and XENON 
indicate a conflict between the experimental results in this case. Alternative couplings, modified 
halo models, and systematic issues have been proposed to reconcile this apparent incompatibility. 
Figure courtesy of J. Kopp |173j . 

impact the interpretation of experimental results. 11 

The high-mass (9(80 GeV) DAMA region appears to be ruled out for both SI and SD 
elastic scattering by null results from CDMS, XENON, and COUPP — a heavy mass WIMP 
is only viable for non-standard interactions. On the other hand, the compatibility of light 
(9(10 GeV) dark matter remains the subject of some debate. For the case of SI scattering, 
these positive results are in apparent contradiction with each other and with CDMS and 
XENON. Some authors have nevertheless argued that some of the results could potentially 
be reconciled (see e.g. Ref. |181j ). 

For the case of SD scattering, the DAMA lower mass region remains compatible with all 
experiments [631 H82[ 1183] , provided the SD coupling to the neutron is strongly suppressed 
relative to the proton (\a a \ <C |a p |). This particular case is uniquely suited to be probed by 
indirect searches involving detection of neutrinos produced by WIMPs annihilating in the 
Sun, e.g. with the Super-Kamiokande |184] and IceCube detectors [185J. 



Considerable discussion remains as to the true sensitivity of the XENON experiment near the energy 
threshold (see e.g. [1751 USD] ). 
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C. Future Prospects 

Direct detection experiments are poised at an important juncture. In the past few years, 
the cross-sections reached by the detectors have improved by roughly two orders of mag- 
nitude. A similar improvement is expected in the next generation of detectors, which will 
be one tonne (1000 kg) in size. These experiments will probe some of the most promising 
regions of WIMP parameter space, exploring Higgs exchange cross-sections and large regions 
of supersymmetric parameter space. However, as the sensitivity of direct detection experi- 
ments reaches <7 Pi si ~ 10 -47 cm 2 , astrophysical neutrinos become an irreducible background, 
so the experiments are no longer zero-background |186H188] . In addition to tonne-size de- 
tectors pushing the reach to lower cross sections and heavier dark matter masses, efforts 
are also being made to explore dark matter with masses below ~1 GeV using electron 
recoils [miWOj. 

New technology and creative experimental designs will allow for further exploration of 
the O(10 GeV) dark matter anomalies. For example, DM-Ice [32] is a detector located 
at the South Pole that uses the same target material as DAMA. Because it is located in 
the southern hemisphere and is embedded deep in the ice where the natural temperature 
variation is minimal, DM-Ice should have different environmental background sources than 
DAMA. DM-Ice should thus be able to conclusively confirm or refute the annual modulation 
claim by DAMA. 

In addition, directional detectors will provide a powerful probe in mapping out the distri- 
bution of the local dark matter. These detectors rely on the fact that Earth's motion with 
respect to the Galactic rest frame also produces a diurnal modulation of the event rate |191] . 
The WIMP flux in the lab frame is peaked in the direction of the Sun's motion and, as a 
result, the recoil spectrum is peaked in the opposite direction for most energies. Therefore, 
the event rate should experience a strong forward-backward ('head-tail') asymmetry. In 
addition, the direction of the dark matter as observed in the lab frame changes with the 
time of day due to the Earth's daily rotation. The modulation amplitude of this diurnal 
modulation can be as large as ~100% |1914 H92J. Ref. [193] reviews the current status of 
prototypes of directional detection experiments. To achieve reasonable angular resolution, 
the recoiling nucleus must leave a track that is sufficiently long. As a result, the chosen 
detector material is a gas, typically CF 4 and CS2 in current designs. The use of gas as the 
active target, with the gas being at low pressure (well below atmospheric pressure to allow 
for longer recoil tracks), will require these detectors to have volumes of O(10 4 m 3 ) to achieve 
tonne-scale masses. 

A novel type of directional detector has also recently been proposed that uses a DNA 
tracking material |194j . These detectors can achieve nanometer resolution with an energy 
threshold of 0.5 keV and can operate at room temperature. When a WIMP from the Galactic 
halo elastically scatters off of a nucleus in the detector, the recoiling nucleus then traverses 
thousands of strings of single stranded DNA (ssDNA) and severs those ssDNA strings it 
hits. The location of the break can be identified by amplifying and identifying the segments 
of cut ssDNA using techniques well known to biologists. Thus, the path of the recoiling 
nucleus can be tracked to nanometer accuracy. By leveraging advances in molecular biology, 
the goal is to achieve about 1,000-fold better spatial resolution than in conventional WIMP 
detectors at a reasonable cost. 

Directional detectors are particularly useful in mapping out the local dark matter distri- 
bution |195H200j . A positive signal at both a direct and directional detection experiment 
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would provide complementary information about the halo, building our understanding of 
the velocity structure of the local dark matter. 

V. SUMMARY 

The theoretical and experimental status of the annual modulation of a dark matter signal 
(due to Earth's rotation around the Sun) in direct detection experiments has been reviewed 
here. Annual modulation provides an important method of discriminating a signal from 
most backgrounds, which do not experience such a yearly variation. The Milky Way halo 
consists of a dominant smooth component as well as substructures such as streams, tidal 
debris, and/or a dark disk, each of which contributes to the modulation of the signal. In the 
Standard Halo Model, the count rate in experiments should peak in June with a minimum 
in December; substructure may change the phase, shape, and amplitude of the modulation. 
The current experimental situation is puzzling, as several experiments have positive signals 
(DAMA and CoGeNT both see annual modulation and CRESST-II has unexplained events) 
but appear to be contradicted by null results from other experiments (CDMS and XENON). 
In the future, detectors with sensitivity to the directionality of WIMPs should enable deter- 
mination of the direction of the WIMP wind as well as diurnal modulation due to Earth's 
rotation. Proposed techniques for directional detection include large gaseous detectors as 
well as nanometer tracking with DNA. Consistent measurement of a head/tail asymmetry 
together with annual modulation would provide very convincing evidence of WIMP detec- 
tion. The future of dark matter searches is promising and an annual modulation signal 
should play an important role in the interpretation and confirmation of a potential WIMP 
signal. 
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Appendix A: Quenching Factor 

Any experimental apparatus does not directly measure the recoil energy of scattering 
events. The recoiling nucleus (or recoiling electron, in the case of some backgrounds) will 
transfer its energy to either electrons, which may be observed as e.g. ionization or scintillation 
in the detector, or to other nuclei, producing phonons and heat; these are the signatures 
that are measured. Some experiments that measure only scintillation or ionization give their 
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results in terms of the electron-equivalent energy E ee of an event in their detector (usually 
given in units of keVee). This quantity is defined as the energy of an electron recoil that 
would produce the observed amount of scintillation or ionization, even if the event was 
actually a nuclear recoil rather than an electron recoil. Nuclear recoils tend to produce a 
smaller amount of scintillation/ionization than electron recoils for the same recoil energy, so 
E ee is not the recoil energy of that event if it is a nuclear recoil event. For nuclear recoils, 
these two energies are related by E ee = QE nr , where Q is called the quenching factor. The 
quenching factor is different for each element in a detector and can have a recoil-energy 
dependence. The different quenching factors for different elements and for electron-recoil 
events (Q = 1 for electron recoils, by definition) makes it impossible to determine the recoil 
energy of an event based upon the scintillation or ionization signal alone. 

Take, for example, the Nal in DAMA, with Q Na « 0.3 and Qi « 0.09 [201] (see Ref. |2D2] 
and references therein for quenching factor measurements of Nal and several other scintilla- 
tors used in direct detection experiments). A recoil event in DAMA that produces 2 keVee 
of scintillation can be from a ~7 keV Na recoil, a ~22 keV I recoil, or a 2 keV electron recoil. 
DAMA is unable to distinguish between these three types of events on an event-by-event 
basis, so any DAMA analysis is necessarily based upon the total 4| spectrum that contains 
contributions from all three types of events. 

CoGeNT, which observes only ionization in a germanium target, also gives results in 
terms of the electron-equivalent energy spectrum. CoGeNT suggests using Q(E xa ) = 
0.19935-E^ 1 as a reasonable approximation to the quenching factor measurements over 
the energy range of interest |203j (measurements of Q in germanium can be found in the 
Appendix of Ref. |47] ) . CDMS, which also uses a germanium target, can discriminate be- 
tween electron-recoil and nuclear-recoil events and can reconstruct the nuclear recoil (though 
limited by a finite resolution), so results for this experiment are given in terms of E ~ E m 
rather than E ee . Because CDMS and CoGeNT are made of the same target material, one 
would expect that these two experiments should have the same nuclear recoil spectrum; how- 
ever, one should keep in mind the caveat that the two results are given in terms of different 
quantities that must be rescaled to make direct comparisons between their two results. 

Appendix B: Mean Inverse Speeds of Commonly Used Velocity Distributions 

The detection rate in dark matter experiments is directly proportional to the mean inverse 
speed rj(v m - m ) = J^ >v d 3 v . Here, we present analytical results of this integration 

quantity for several commonly used isotropic velocity distributions. We define /(v) as the 
velocity distribution in the rest frame of the dark matter population (i.e. J d 3 uv/(v) = 0). 
The velocity distribution in the lab frame is determined via the Galilean transformation 
/(v) = /(v bs + v), where v Q b s (t) is the (time-dependent) motion of the lab (observer) 
relative to the rest frame of the dark matter population. This motion is described in more 
detail in Section IHII 

The distributions considered below are Maxwellian distributions (including two modifica- 
tions to account for a finite cutoff) and distributions corresponding to cold flows and debris 
flows. These are not the only possible distributions of dark matter and may be only simple 
approximations for some populations, but they are frequently used distributions that have 
known analytical forms for rj. 
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1. Maxwellian Distributions 



Perhaps the most useful simple distribution is the Maxwellian: 
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where v is the most probable speed, 

x = v min /v and y = v ohs /v . 
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Many well-mixed populations of dark matter particles can be expected to have a Maxwellian 
or Maxwellian-like distribution, in which case the above is a useful first approximation. 

The Standard Halo Model (SHM) takes the dark matter halo to be an isothermal sphere, 
in which case the velocity distribution is Maxwellian. However, high velocity particles would 
escape the Galaxy, so the high-velocity tail of the distribution is cut off in a realistic halo 



model. Section II B 1 presents two methods for removing the tail of the Maxwellian, with 



the resulting distributions given by Eqns. (14) & (17). For the SHM or any dark matter 



component described by one of these two velocity distributions, the mean inverse speed rj is 
[T2T1 [155] 
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where /3 — for Eqn. (14) and /3 — 1 for Eqn. (17), x and y are as defined above, and 



z = v esc /v . 

Note that the normalization factor N P v C has a different form for the two distributions. 



(B4) 
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2. Cold Flow 

Cold flows, such as tidal streams, have small to negligible velocity dispersions. In the 
case of zero dispersion, 

/(v) = 5 3 (v) (B6) 
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and 

?7(Vmin) = #(^obs - *W) , (B7) 

^obs 

where 9 is the Heaviside function. Although the velocity dispersion is small in cold flows, 
in some cases such as tidal streams, it is not completely negligible. In those cases, a useful 



approximation can often be made with the Maxwellian distribution of Eqn. (Bl) with a 
small Vq. 



3. Debris Flow 

To first order, the debris flow in a Milky Way-like galaxy [126] has an isotropic, constant- 
speed velocity distribution in the Galactic rest frame that can be described by a delta 
function in speed, as opposed to the delta function in velocity seen with cold flows above: 



/>) = ^^(|v|-tw), (B8) 



flow 

where t>fl ow is the uniform speed of the particles. In the lab frame |123j . 
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Although the debris flow has both a dispersion in the speed and some anisotropy in the 
Galactic rest frame, the above reference shows that this analytical form still provides a 
reasonable approximation. 
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